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ABSTRACT 


The algorithms for processing CZCS data to geophysical units (pigment 
concentration) are described. This document provides a summary of current 
public-domain information for processing these data. Calibration, atmospheric 
correction, and bio-optical algorithms are presented. Three CZCS data- 
processing implementations are compared. 
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DESCRIPTION OF ALGORITHMS FOR PROCESSING CZCS DATA 
I. INTRODUCTION 

■r* 

— i 

This document Is a short Introduction to the Coastal Zone Color Scanner 
(CZCS) data and a description of the extant versions of the CZCS algorithm. 

The CZCS is a scanning multispectral radiometer with a scan angle of 78 
degrees and an instantaneous field of view (IFOV) of 0.0495 degrees. The 
nominal resolution at 955-kra altitude is 825 meters at nadir. To avoid sun 
glint, the sensor can be tilted forward or backward in 2-degree increments, to 
a maximum of 20 degrees. The four bands on the CZCS relevant for color work 
are 20 nm-wide bands centered at 443, 520, 550, and 670 nm. Another band at 
700-800 nm is used for detecting land or cloud pixels. A sixth band in the 
emitted infrared has functioned intermittently during the mission. 

The CZCS scans 1968 samples per line. There are two major formats for 
digital data on magnetic tape, the ZIP format (e.g., Scripps or RSMAS tapes) 
and NESDIS Level-1 tapes.* This document does not address tape formats or 
unpacking of data. 

The "color algorithm" is threefold (see Figure 1). First, the sensor- 
apparent total radiances, L T , are derived from the satellite digital counts 
using a calibration algorithm. Then the water upwelling radiances, L^, are 
derived from the values using an atmospheric correction algorithm. 

Finally the water upwelling radiances are used to derive the pigment concen- 
tration using a "bio-optical" algorithm. 


*Scripps is the Scripps Institution of Oceanography Satellite Facility. 
RSMAS is the Rosenstiel School for Marine and Atmospheric Sciences at the 
University of Miami. NESDIS is the National Environmental Satellite Data 
and Information Service. All three have various degrees of CZCS data 
processing capability. 



Figure 1. Generic CZCS Algorithm (after Gordon et al., 1983a) 










II . CALIBRATION ALGORITHM 

During data acquisition, sensor-apparent radiances are digitized for 
transmission to earth* Calibration equations are used to transform satellite 
digital counts back to sensor-apparent radiances. Following are the cali- 
bration procedures and constants necessary for the conversion. Also included 
are some other constants needed for later processing. 

Sets of calibration coefficients are chosen according to the gain setting 
of the CZCS (possible range 1-4) for a particular pass. The gain setting must 
be available to any calibration program. 

The equation to convert counts to radiances can be of the form: 

(A(A)N(A) + B(A)]C'(A) = L^( A) (the total radiance) 

where 

N * digital count from the sensor 

A = calibration slope 

B = calibration intercept 

C = sensor response correction function, f(C). Determination of C 
depends on the value of C, a correction coefficient. C in turn 
depends on estimates of the mean extraterrestrial solar radiance 
F n . In short, these three parameters are coupled. Care must be 

taken to use coherent sets of C', C, and F^. See e.g. Viollier 
(1982) for other coefficients. The values given here are from Austin 
(quoted in Gordon et al. (1983a)) and have the advantage of being 
independent of gain. 


A * wavelength of one of the four visible bands sensed by the CZCS. 



It is important to note that the CZCS sensor response has been changing 
with time (Gordon et al. (1983b)). Current efforts in CZCS calibration 
therefore define C'(A) as a function of orbit number (time). This sensor 
degradation correction is the focus of much current activity in the CZCS 
community. 

The coefficients for CZCS parameters A and B in Table 1 were furnished by 
Mueller (personal communication). The figures for F q , m, and the atmospheric 
optical parameters in Table 2 are from the Goddard Space Flight Center (GSFC) 
(undated, but probably 1982). All of these are generally undisputed. Figures 
not universally accepted at the time of this writing are given with references. 
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Table 1. CZCS Coefficients 


L 




f 

£ 

■h; 

443 nm 

520 nm 

550 nm 

670 nm { 

• 


\ (in mW/cm^-vtm-sr-count ) 


Gain 1 

0.04452 

0.03103 

0.02467 

0.01136 

Gain 2 

0.03589 

0.02493 

0.02015 

0.00897 

Gain 3 

0.02968 

0.02032 

0.01643 

0.00741 

Gain 4 

0.02113 

0.01486 

0.01181 

0.00535 

i mW/cm^- 

■|im-sr ) 




Gain 1 

0.03963 

0.06361 

0.07992 

0.01136 

Gain 2 

0.05276 

0.08826 

0.06247 

0.03587 

Gain 3 

0.02879 

0.09752 

0.06570 

0.02963 

Gain 4 

0.03359 

0.05647 

0.04723 

0.01604 

ordon et 

al. (1983a), GSFC) 





1.0688 

0.9931 

0.9554 

1 .000 

ueller ) 

1.144 

1.033 

0.979 

1 .000 


Sensor Degradation Correction Functions, C': 

Gordon et al. (1983b): 

C'(X) = C(X)[a - b(orbit) + c(orbit^)] 
for X * 443 nm, a = 1.086, b = 2.46 x 10 c 

a = 1.069, b = 2.32 x 1(T 5 , c 
(see reference) 

for X ■ 520 nm, a = 1.024, b » 0.59 x 10~\ c 

for X = 550 nm, a = 1.007, b = 0.28 x 10"^, c 

Mueller : 

C'(X) = C(X)exp[a(orbit)] 
for X = 443 nm, a = 2.12 x 1U - ^ 

for X = 520 nm, a = 1.22 x 10“^ 

for X = 550 nm, a = 0.78 x 10 - ^ 


5.05 x 10“ 
5.0 x 10" 1 




C'(A) = C(A)expla(orbit - 3200)] 

for A = 443 nm, a = 2.12 x 10~^ 

for A = 520 nm, a = 1.22 x 10~ 5 

for A = 550 nm, a « 0.78 x 10"^ 

F 0 (in mW/cm^-um) 

186.42 185.34 184.76 


m (nominal refractive index of seawater) 






Table 2. Atmospheric Optical Coefficients 



443 nm 

520 nm 

550 nm 

670 nm 

R ” Rayleigh optical thickness 
< 25° Lat. 0.2329 

0.1231 

0.0969 

0.0444 

2 50 -550 Lat. 

Summer 

0.2311 

0.1222 

0.0962 

0.0440 

Winter 

0.2316 

0. 1224 

0.0964 

0.0442 

> 55° Lat. 

Summer 

0.2300 

0. 1214 

0.0956 

0.0438 

Winter 

0.2303 

0.1218 

0.0959 

0.0439 

oz ~ Ozone optical thickness 
< 25° Lat. 0.0066 

0.0166 

0.0261 

0.0158 

2 50-550 Lat. 

Summer 

0.0067 

0.0200 

0.0323 

0.0191 

Winter 

0.0069 

0.0237 

0.0390 

0.0226 

> 55° Lat. 

Summer 

0.0068 

0.0213 

0.0346 

0.0202 

Winter 

0.0071 

0.0275 

0.0461 

0.0264 


Notes : 


(1) Atmospheric optical parameters are dependent on the time of 
year and the latitude. 

(2) Summer and winter are defined in terms of the local season, 
the summer /winter boundaries being the spring and autumn 
equinoxes . 
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III. ATMOSPHERIC CORRECTION ALGORITHM 


» . 





The following description is of a color algorithm taken from the Ocean 
Color Science Working Group (OCSWG) (1982) and Gordon et al. (1983a). Further 
references for the interested reader may be found in those papers. The detail- 
ed provenance of this algorithm is not known; the term "Gordon algorithm" is 
used as a convention in this document because he is first author of the 1983 
paper that sets out clearly in the open literature the steps necessary to 
implement a generic color algorithm and because he first suggested the general 
approach (Gordon, 1978). Many of the equations here are taken from Gordon et 
al. (1983a) and Sturm (1981). Other algorithms do exist, however. The Gordon 
algorithm is the basis of the GSFC production processing, the RSMAS procedures, 
and the Mueller procedures. Differences among the RSMAS, Mueller, and GSFC 
procedures are discussed in Section V, Implementations. 

The total radiance apparent to the CZCS may be partitioned into conti i- 
butions from sun glint (specular reflection from the sea surface), Rayleigh 
and aerosol backscattering, and the water-leaving radiance. The latter carries 
information about the pigment concentration of the water. The CZCS mirror is 
tiltable to avoid glint. The first step in the color algorithm is removing 
the atmospheric contribution at each wavelength from L^. Assuming no glint, 

Lj is given by: 

h T (X) = L r (X) + L a (X) + t(X)L w (X) (1) 


where 

L * radiance due to Rayleigh scattering 
R 

= radiance due to aerosol scattering 
= radiance leaving the water 
t * diffuse transmittance to top of the atmosphere 

Since calculations of pigment depend on L»j, we shall solve for that term. 

We know from physics by way of Gordon et al. (1983a) that: 

u .a) T A U)F '(A)p.(e,e .x) 



a) = single scattering albedo of the aerosol (good approximation 

A 

to multiple scattering) 
t a = optical thickness of the aerosol 
p, v (6,0 o ,X) = aerosol characterization function, considering the Fresnel 

reflectance of the interface and the aerosol scattering phase 


function 

F ' = spacecraft apparent Instantaneous extraterrestrial radiance, 
o 

F , corrected for two passes through the ozone layer, a known 
o 

constant. In practice, F 0 values are the mean extra- 
terrestrial radiation values, F () , corrected for the CZCS 
sensor response and the time of year. These values must 
be highly consistent with the sensor calibrations listed in 
Section II, Calibration Algorithm, ar.d they have been so refined 
(see Gordon (1981) for discussion of F 0 determination). 

The F ' values are related to F by: 
o ° 



[1 + 0.0167 cos ( — 3 ^ 5 ~ ) 1 ^ 
exp 1 " T 0Z ^cos“e + cosT^ 1 


where 

D = Julian day of the year ( 1< D < 366) 

9 = solar zenith from pixel local tangent plane 

9 = spacecraft zenith from pixel local tangent plane. Zenith angle 
measurements must be based on knowledge of sun, spacecraft, and 
pixel location. Therefore the image must be navigated, i.e., 
the pixels must be earth-located. (This also applies for azimuth 
angles that appear later in the algorithm.) A discussion of 
navigation is beyond the scope of this document. However, in 
general, navigation depends on knowledge of the orbital charac- 
teristics of the spacecraft, accurate time of overpass, and 
knowledge of spacecraft attitude during overpass. There may be 
some ambiguity in location description. A specific example of 
configuration is given in Figures 2 and 3. This configuration 
is used for Equations 7—10. General illustrations are given in 
Section V.C, Figures 5-10. 
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If we assume a constant aerosol "type", i.e., an atmosphere In which the 
normalized size frequency distribution and the refractive index of the aerosol 
are independent of horizontal position in the image, then for two wavelengths 
A and A^ f we can define S(\ 9 \^) 9 the ratio of L^(A) to l a^o^ # 
ratio S(A > A o ) will be independent of horizontal position in the image, 
even though L A (A) and L a (A q ) may vary. 


If we expand the L.(A) terms in the ratio S(A,A ) we have: 

o 

0) (A)t (A)F '(X)? .(0,6 ,A) 

1 ’ * A ( X )t(a“)f^a )P A (e,e , a ) 

AoAoo oA oo 


(3) 


We can name e (A,A ) as the ratio of the products of the aerosol 
albedo, optical thickness and aerosol function at the respective wavelengths. 
Of the three unknowns, Gordon et al. (1983a) assume that for aerosol, the 
albedo and phase function are weak functions of wavelength or viewing angle, 
and therefore e(A,A^) is determined mostly by the ratios of the 
aerosol optical thickness at the two wavelengths in question. For a single 
aerosol type this ratio should be fairly constant, even if the aerosol is not 
horizontally homogeneous. Then we can state: 


l a< a > 

VV 


S(A,A q ) = 




o 0 
F '(A ) 


where e(A,A^) is assumed constant over the image. 


(4) 


Now substituting into Equation 1, we can solve for the water-leaving 
radiance L W (A): 


tU^A) = L t (A) - L r (A) (5) 

- SU,A o ) |L t (» 0 ) - L r (A o ) - t(» 0 )I v (X 0 )l 

The wavelength A q is chosen so that L^(A q ) is minimized. In this 
way one of the unknowns in Equation 3 is eliminated. In practice, A q is 
the 670 nm CZCS band. Let us consider the terms of this equation left to 
rLght. The diffuse transmittance t(A) is given by: 
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( 6 ) 


-Hr + t > 

t(A) = t A ( A ) exp [ ° z 

A COS 0 

where 

X R ~ Rayleigh optical thickness, a known constant 
T oz = ozone optical thickness, a known constant 
t A = a function depending on the aerosol optical thickness 

In practice, t r and t qz are known. Gordon et al. (1983a) suggest setting 
t A (A) to unity since certain assumptions in the algorithm will break down 
before the aerosol optical thickness starts to affect the diffuse transmittance 

The next known term in Equation 5 is L T , the total sensed radiance. 

This is derived from the digitized counts by certain calibration equations 
(see Section II). 

The next known term is the radiance due to Rayleigh scattering, r .r 
is given by: R 


V A > 


3t oz <A,, » (>,F - <A) 


R \ * \ / 

° / r , 2 0 

( [ 1 + cos Y ] + [ 1 + 


1 6n cos 6 


_] + [!+ cos <f + ] [p(0) + p(0 ))) 


where 

cos i|>_ = [-cos 6 cos 0 Q - Sin 6 sin 0 q cos ($ - $ - n)] 

cos V = [cos 0 cos 0 + sin 0 sin 0 cos (<b -6)1 

O O r O r 7 1 

p(x) = Fresnel transmittance 


(7) 

( 8 ) 
(9) 


= 1 - [ 2m( A )yz ] 


cos x 


(x = 0 or 0 ) 


( 10 ) 



oz^ “ ozone transmittance to top of atmosphere 


= exp -t (A) [ — + — ! j 

oz cos B cos 0 J 


<f> o “ solar azimuth 


* -- spacecraft azimuth angle from pixel to sensor (see Figures ? and 3) 


where 

mv A) - the refractive index of seawater at A 

= tra(X) 2 + cos 2 x - 11 
m( A ) 

z = [cos x + ym(A)]“ 2 + [y + m(A) cos x]“ 2 

These equations have been cast in the form of Sturm (1981). The reader 
should refer to Figures 2 and 3 for proper orientation, as differences in 
angle-naming convention have been noted in the literature. 

The next term in the equation is S(A,A q ) which we shall put aside 
for the moment. 

The next two terms in the equation are derived for A in the sane 
manner as described above for A. Rationale for determination of A - 
670 nm is discussed in the next paragraph. 

The last term is t(A o )L w (A o >. In order to proceed, the atmo- 
spheric correction depends on the quick attenuation in the ocean of radiation 
at 670 nm (CZCS band 4). Assuming no contaminating source in the water column, 
such as sediment from the rivers or phytoplankton, the ocean water itself 

should radiate little at this wavelength. Therefore for clear water L lt (670) = 0. 

w 

(Smith and Wilson (1981) provide an iterative approach to this problem for 
L w (670) * 0.) So if we select 670 nm as A , t(A )L.,(A ) drops out, l.e., 

all the radiance in this band is due to atmospheric effect (for clear water), 
and we can solve for S(A, A ). 


( 11 ) 

( 12 ) 



So, for CZCS wavelengths 443, 520 and 550 nm: 


L (443) - L (443) - S(443,670) [L T (670) - L R (670)] 

V 443) OW 

l. ( 520) - 1^(520) - S(520, 670) [1^(670) - 1^(670)] 
L W (520) * t(520) 

L (550) - L (550) - S(550,670) [L T (670) - L R (670)] 
L w (550) 1(550) 


(i: 


(i< 

(i 


Note the only unknown term on the right side of Equations 13-15 is 
S(A,670). To solve for S(A,670) we must select a clear water pixel. These 
pixels may be assumed to have a known radiance. For clear water pixels with 
pigment concentrations less than 0.25 mg/m , Gordon and Clark (1981) have 
shown that: 


§ + x 

2 02 


L w (A) = [l w (X)1 n cos e o exp 1 ^os“e — 1 


where 


[L v; (520)] n = 0.498 tnW/cm -um-sr, the normalized water-leaving 
radiance at 520 nm 


[ L^ J ( 550 ) ) N = 0.30 mW/ cm" -um-sr, the normalized water-leaving 
radiance at 550 nm 


[ L VJ ( 670) ] N = <0.015 mW/cm"-iim-sr, the normalized water-leaving 


radiance at 670 nm 


(1 


0 


Thus all the terms except S(A,670) in Equations 14 and 15 are known and the 
expressions can be solved algebraically for e(A,670). e(443,670) is 

solved by extrapolation as it has been shown that e(A,670) is a smooth 
function of n (see below). One solves according to the following equations 

e(443,670) - [443/670 ) n(443 ) ( 


o 


n(443) = 0.5 [ 108 c (520.670) log e (550.670) 

1 log (520/670) + log (550/670) 1 

In practice, either a skilled user will select (based on oceanographic 
experience) a clear water pixel from a preview image or an automated procedure 
can be used. One possible automated method is to use a first iteration of e 
to derive trial pigment concentrations. In general, the pixel with the lowest 
pigment value (subject to certain other tests) is acceptable as the clear 
water pixel. Gordon et al. (1983a) spell out the procedure as follows. 
Determination of e's is effected by carrying out the above process using, as 
a first iteration, e(A.) = 1 for i = 1-3. The region of lowest C 
(pigment concentration) is located using Equation 18, then Equations 4 
and 13-17 are solved for e^, e^, and e^ using values averaged 
over a 5 x 5 pixel region to minimize ncise. The "region of lowest C" 
must also meet the test of having a monotonically decreasing sequence 

La 52o’ La 550* 3nd La 670* ^ theSe ’ the region with the highest is 

selected. This region is used as the clear water calibration area, e's are 
derived using Equations 13-17 and processing continues using the bio-optical 
algorithm throughout the scene. 


In summary, the use of clear water pixels as areas of known Ly allows 
the calculation of the atmospheric parameter e, so that correct L„ value-: 

W 

can be computed throughout the scene without simultaneous in-water radiance 


measurements. 


Once a satisfactory set of e's is derived, these are applied for each 
pixel in the scene, and the usual route of extracting L^A) from L r[ .(A) 
is followed throughout the image, by Equations 13-15, using S(A,A ) 
derived from e(A,A ). 


y 


Gordon et al. (1983a) give procedures for atmospheric correction in three 
cases that violate the assumption of a constant aerosol type. The three cases 
are: variable aerosol over clear water, two distinct aerosol types with the 
more turbid type overlying at least some clear water, and two distinct aerosol 
types with the more turbid type not overlying some clear water. The first 
case is remedied by deriving n(A) at each pixel. The second is mechanistically 
identical to the usual procedure as described above. The third case can be 
handled by manually varying values of n(A). This requires somewhat different 
procedures to accomplish atmospheric correction in at least a semiautomated 
fashion. 



IV. BIO-OPTICAL ALGORITHM 


The bio-optlcal algorithm is a straightforward transformation from water 
upwelling radiances, Ly, to pigment values (see Clark, 1981). ("Pigment" 
neans chlorophyll and its associated degradation products. The spectral 

signatures of chlorophyll and its degradation products are indistinguishable 
by the CZCS.) 


Sea-truth data measurements have yielded the following regression 
equations for conversion from water upwelling radiances to pigment 
concentration in mg/m 3 . 


log 10 C = 0.053 + 1.71 


Vi 


1°810 


550 


W 


443 


(is; 


log io C = 0.522 + 2.44 


logio 




(19; 


Equation 18 is used as the default. If C derived from Equation 18 is more 
than 1.5 mg/ny then C is determined by Equation 19 if it produces C greater 
than 1.5 mg/m . This is because when C is high, absorbance by pigment 

causes L to become too small to be retrieved with sufficient accuracy 
W 443 

from Lj. . The alternative Equation 19 does not use L 

443 W, ... 


These coefficients have been tuned to L w , the water-leaving radiance. 
The GSFC procedures convert water-leaving radiance (i.e. just above the 
air-sea interface) to subsurface radiance, L gg (i. e . just below the air-sea 
interface), by correcting for Fresnel reflectivity losses through the 

interface (L gg » m Ly(A)/l - p(0)). That algorithm uses different 

coefficients in Equations 18 and 19, however. According to Gordon (personal 

communication), there is no significant difference in using the ratios of the 
U.'s vs the L 's. 

w S3 


V. 


IMPLEMENTATIONS 


A. RSMAS 

The University of Miami's Rosenstiel School of Marine and Atmospheric 
Sciences (RSMAS) has implemented a software system for processing CZCS data. 
General features of this system are described in Brown and Evans (1982). 
According to Gordon (personal communication) the RSMAS implementation is as 
described in Gordon et al. (1983a). RSMAS is conducting research on changes 
in CZCS sensor response, among other things. 


B * Dr. J. Mueller, Naval Postgraduate School 

Figure 4 is a flow diagram of the Mueller algorithm. Dr. Mueller has 
given us code for processing a navigated image. References to subroutines are 
those listed in Figure 4. It may be instructive to consider the modular nature 
of Mueller's implementation to gain insight into the required processing steps. 

In general, this algorithm does not calculate each parameter at each pixel 
but rather at every 16th pixel ("anchor point" in CZCS terminology). This 
approach is also mentioned in Gordon et al. (1983a). To provide the 

parameters for each pixel, gradients are calculated between the anchor points 
in the x and y directions. 


The program MAIN calls some routines particular to the NPS environment. 
These routines are JOBNAM, FILNAM, RUNLOG, and TOD - they need not concern us 
Another subroutine called by MAIN is GLOAD. GLOAD gets the navigation para- 
meters, sun zenith and azimuth, and spacecraft zenith and azimuth, for the 
anchor points. These parameters are copied to a direct-access file for later 


MAIN Chen cell. SETUP. SETUP eete up some erreys and common blocks. It 
reads in a NAMELIST of parameters for defining a data window, reads the gain, 
reads the exponent for Equation 17, and reads a threshold for land/cloud flag, 
One difference between the Mueller and Gordon Implementations Is that Mueller 
depends on an estimate of n(A) (see Equation 17) to derive t(X), while 
Gordon uses estimates of . derived from the aerosol radiances to estimate 
n(A). This is a conventional rather than a substantive difference. Other 
parameters sec flags for printing some diagnostic output. SETUP calls sub- 
routines ABEND and RUNLOG, both site-specific routines. It also calls three 
routines, SUNFLX, SETCAL, and OPTICS to read In more coefficients. SUNPLX 
calculates the F o values (see Equation 2 for F o ). SETCAL reads in 

calibration coefficients. OPTICS reads in some atmospheric optical 
coefficients . 


The next routine called by MAIN is PROCES, which does line-by-line proce 
smg (i.e., it is counting lines). PROCES in turn calls SETBLK (sets up the 
arrays for interpolation), GETLIN (gets a scanline of data from tape), LINCAL 
computes radiances and flags land/cloud pixels), and WRTLIN (outputs a line 
of data, either radiance data or data flagged as bad). 

SETBLK calls GLINES and GREAD, which retrieve a navigation record for the 
current line from a direct-access file. SETBLK then calls, through ORGVAL 
GEOMTY and RADCAL. GEOMTY calculates the angles and Y + and atmo- 
spheric slant paths for later use. RADCAL calculates L R ' S and some inter- 
mediate parameters for calculating L w 's. SETLLK then calls GRADS, which 

through XGRAD and YGRAD, fills the interpolation buffers and repeats the GEOMT 
and RADCAL calculations for each pixel. 

GETLIN reads in a line of data and does some input checking. 

LINCAL does the conversion from counts to radiances. First it flags 
Pixels as land/cloud if the band 5 count is greater than a predetermined 

The remaining pixels are processed in subroutines PIXCAL, PIXOUT 
and CLPIX. PIXCAL derives I^'s and then pigment and k, the diffuse attenua- 
tion coefficient. PIXOUT fills a data buffer, including the unprocessed band 
counts. CLPIX is an alternative output routine for diagnostic purposes. 


The fourth subroutine called by PROCES is WRTLIN, which does the physical 
output of data. 

The remaining subroutines in Figure 4 are DIAGNS and FCON. DIAGNS is a 
routine to dump some diagnostic output if a certain flag is set. FCON is a 
common block used by many routines. 

A few unique features of the Mueller algorithm are its use of bilinear 
interpolation, the estimate of n(A) rather than e in Equation 17, and the 
processing of data on a swath basis rather than on subsets of a swath. 
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c. 


GSFC Production Algorithm 


The GSFC production algorithm is used by NASA to generate products for 
use by the public. These products are available through NESDIS. They are the 
closest thing to operational products that exist for CZCS. 

The GSFC algorithm depends on partially preprocessed data and ancillary 
information. This information is contained in the so-called Level-1 tape, 
which is the input for the GSFC procedures. This information is needed for 
navigation of the image. Consideration of which data are used by the GSFC 
algorithm may be instructive. 

The first step in the GSFC algorithm is calibration of counts to radi- 
ances, including the sensor degradation correction function. This has been 
presented in Section 2. 

The next step is the geometry calculation. This takes place in four 
stages: scene global, scan line, anchor point, and final geometry calcu- 

lations. The geometry calculation depends heavily on vector arithmetic for 
the navigation. 

The scene global calculation takes place once per scene. The range over 
the scene of ephemeris time, sun position, spacecraft position, and Greenwich 
hour angle are calculated. Ephemeris times are on the Level-1 tape. The sun 

vector is determined from sun right ascension and declination data on the 
Level-1 tape. 

The scan line stage calculates for each line the ephemeris time, sun 
position, spacecraft position, and hour angle. These calculations are just 
interpolations in time of the global parameters calculated above. 

The anchor point stage uses some 77 predefined anchor points. The Level-1 
tape contains preprocessed latitude/longitude coordinates for these locations. 
The coordinates are converted to vectors. This stage then computes a 
spacecraf t-to-pixel unit vector for each point. 



0J 


The final geometry calculation computes for each anchor point the solar 
zenith angle 0^, the spacecraft zenith angle 0, the scattering phase 
angles (for use in deriving Rayleigh radiance), the solar azimuth angle <j> o , 
and the spacecraft azimuth angle <\> (see Figures 5 through 10)* These values 
can then be linearly interpolated in the direction of scan for the non-anchor 
point pixels, providing a full set of geometry values for each pixel. 

After navigation the procedure is similar to the Gordon algorithm (see 

Figure 1). Solar and atmospheric constants are derived from NAMELISTs. To 

derive e(A)'s, the GSFC algorithm searches alternate pixels on alternate 

scanlines to locate a likely clear-water area. There are three criteria used 

to test for acceptable clear-water areas. They are spacecraft and solar zenith 

2 

angles >0.6 radians, radiance at 670 nm less than 1.4 mW/cm -sr-um, and 

\j /[L. + L ] between 0.9 and 2.0. Then calculation of the clear- 

T 443 T 520 T 550 

water radiance can proceed according to Equations 16-16c and e(A) can be 
derived by Equations 13-15. 

One difference between the GSFC and Gordon algorithms is the two coeffic- 
ients in Equations 16a and 16b. The GSFC algorithm uses 0.495 and 0.280 
respectively, rather than Gordon's 0.498 and 0.300. Another slight difference 
is in the method of calculating Rayleigh radiance. The two methods also differ 
in selection of a clear-water area. The GSFC algorithm selection examines the 
calculated sets of e values. Sets of e(A,A^), with a sequence of 
e(520,670), e(550,670), 1 which is not monotonic are excluded as are set** 
with £(443, 670) greater than 3.0. Unique to GSFC is the method of choosing 
e( A, A q). The remaining sets of e(A,A Q ) are compared, and the 
set having the lowest e(443, 670)/L,p(670) is chosen as the best for atmo- 
spheric correction. 
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Figure 6. Solar Zenith Angie (from GSFC, 1982) 






Figure 9. Solar Azimuth Angle (from GSFC, 1982) 



Figure 10. Spacecraft Azimuth Angle (from GSFC, 1982) 
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